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Введение. Одной из особенностей обработки трубных досок большого диаметра на радиально- 
сверлильных спецстанках для сверления глубоких отверстий на основе станков типа «Чепель» 
является возникновение в процессе обработки автоколебаний системы СПИД в вертикальном на- 
правлении. Возникновение автоколебаний является следствием относительно малой жесткости 
как преобразующей системы станка, так и обрабатываемой заготовки [1]. Практика эксплуатации 
описанного станочного оборудования показала, что возникающие автоколебания, в зависимости 
от частоты и амплитуды, могут иметь и благотворное влияние на процесс обработки и негативное. 

Цель исследования — повышение технологической надежности процесса сверления отвер- 
стий (в том числе, глубоких) в трубных досках большого диаметра из труднообрабатываемых ма- 
териалов, а также в трубных досках из слоистых материалов на спецстанках, изготовленных на 
основе радиально-сверлильных станков «Чепель» с использованием твердосплавных сверл с 
внутренней подачей масляной СОЖ. 

Поставленная цель может быть достигнута в результате решения комплекса объемных на- 
учных и технологических задач, составляющих содержание исследования: 

— разработка динамической математической модели подсистемы инструмента спецстанка; 

— разработка динамической математической модели подсистемы заготовки и рабочего 
стола спецстанка; 

— разработка математической модели процесса сверления глубоких отверстий, учитываю- 
щей особенности сверления труднообрабатываемых материалов; 

— построение объединенной динамической математической модели, включающей в себя 
вышеупомянутые модели; 

— идентификация параметров математической модели по результатам экспериментальных 
исследований, проведенных при сверлении образцов трубных досок, а также по результатам не- 
посредственных измерений и расчетов; 

— аналитическое и численное, с использованием методов имитационного моделирования, 
исследование динамических свойств объединенной модели. Определение возможности и условий 
возникновения автоколебаний в системе сверления; 

— разработка методов выбора режимов обработки для получения возможности управления 
автоколебательными режимами при глубоком сверлении труднообрабатываемых материалов. 


Работа выполнена при поддержке ФЦП «Научные и научно-педагогические кадры инновационной России» 
(ГК 14.740.11.1323) 
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В ходе исследования определены: 

— возможность и условия возникновения автоколебаний в модели системы сверления; 

— наблюдаемость и управляемость модели по управляющим воздействиям и параметрам 
системы управления; 

— области устойчивости и области возникновения устойчивых автоколебаний в динамиче- 
ской системе сверления. 

В исследовании представлено решение первой из поставленных задач — построение дина- 
мической математической модели станка для сверления глубоких отверстий. 
Математическая модель преобразующей подсистемы спецстанка. Внешний вид радиаль- 
но-сверлильного станка КЕ-51 «Чепель», на базе которого создан спецстанок для сверления глу- 
боких отверстий, приведен на рис. 1. Станок состоит из колонны, траверсы, сверлильной бабки. 
На рабочем столе расположен обрабатываемый пакет перегородок пароконденсатора. Упрощен- 
ная кинематическая схема станка приведена на рис. 2. 
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Рис. 2. Кинематическая схема модели станка 
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Математическая модель, включающая в себя все кинематические подсистемы станка, по- 
лучается достаточно объемной (станочная система имеет не менее 7 степеней свободы). Однако 
для анализа устойчивости системы и определения закономерностей возникновения автоколеба- 
ний, связанных с взаимовлиянием подсистем через процесс обработки, она может быть сущест- 
венно упрощена. В частности, динамическими процессами в зацеплениях между червячной пере- 
дачей и зубчатым колесом можно пренебречь без изменения качественного поведения системы, 
сведя их к эквивалентному моменту инерции .Л/‚,. Также можно пренебречь динамическими про- 


цессами в редукторе между приводом главного движения и шпинделем, рассматривая редуктор 
как идеальное усилительное звено. Подсистему закрепления сверлильной бабки на траверсе 
можно заменить эквивалентной упруго-диссипативной системой, учитывая изменение плеча ры- 
чага переменными коэффициентами. Преобразованная эквивалентная кинематическая схема 
спецстанка показана на рис. 3. 
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Рис. 3. Упрощенные кинематические схемы подсистемы привода 
главного движения (а) и подсистемы подачи (6) 


Подсистема привода главного движения. Механическая подсистема вращения шпинделя 
характеризуется углом поворота Вала двигателя о, и скоростью вращения «,, вращательным 


моментом М,, приложенным к ротору, моментом инерции ротора электродвигателя и упругой 
полумуфты сочленения со шпинделем станка ./,,, коэффициентом диссипации в подшипниках 
двигателя й, ,, жесткостью и диссипацией в упругой муфте с, и й, , углом поворота шпинделя с, 
и скоростью его вращения ‹,, моментом инерции шпинделя станка .Л,, коэффициентом диссипа- 


ции в подшипниках шпинделя /,, и моментом сопротивления вращению, формируемому процес- 


сом резания М, . Система уравнений, описывающая вращение шпинделя, имеет вид: 
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Подсистема сервопривода. Упруго-диссипативная подсистема с коэффициентами жесткости с, 
и диссипации й обеспечивает передачу усилия от сервопривода на шпиндель, обеспечивая его 


поступательное перемещение. Принятые обозначения: й,, — коэффициент диссипации в подшип- 


никах; х -— координата подвижной части сервопривода; у» - скорость; х, — координата поступа- 
тельного перемещения шпинделя; у, — скорость; т, - эквивалентная масса подвижной части 
сервопривода; ти, — эквивалентная масса подвижной части шпинделя; Е - усилие, действующее 
на сервопривод со стороны серводвигателя [2]; Р’ - усилие, действующее на шпиндель со сторо- 


2. 


ны процесса обработки. Влияние траверсы учитывается коэффициентами жесткости и диссипации 
с, И й, и эквивалентной массой т„,. Демпфирующий коэффициент 7,, учитывает потерю энер- 


гии в пластинчатой пружине в подсистеме подвеса шпинделя. 
Уравнения динамики подсистемы имеют вид. 
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Моделирование процесса резания. Процесс резания выступает нелинейным элементом, объе- 

диняющим подсистемы сверлильного станка. Поэтому построение корректной модели процесса 

сверления является принципиально важным для достижения поставленной цели исследования. 
Сила резания в первом приближении может быть описана [3] выражением: 


Е, —- РЕЗ 5 ! 
где р - коэффициент резания; {, — ширина срезаемого слоя; 5„ — оборотная подача. 
Осевая составляющая силы резания в этом случае составляет РЁ, =/..рЁ.5„, момент со- 


противления вращению М, =Л.рё5 „г, где г -— радиус сверла. Учитывая, что для спиральных 


сверл 1, =г, имеем: 
2 
Е =А рожи, М, =А РЭ’. (3) 
При равномерном сверлении 5. определяется отношением скорости вертикальной пода- 


чи к частоте вращения сверла (в единицах СИ): 
У 
> 4 
б Г. 
1800; (9) 
т. е. зависимость оборотной подачи от скорости подачи и частоты вращения формирует нелиней- 
ную связь между подсистемами станка через силовые характеристики процесса резания: 


пу, › ПУ, 
Ао В. Я 
1806, г ^:Р 1806, (5) 











М, =А.р 


65 


Технические науки 


Для понимания процессов, приводящих к заклиниванию при глубоком сверлении, принци- 
пиальное значение имеет взаимовлияние приводов [4, 5]. При этом можно пренебречь подсисте- 
мой траверсы и промежуточных упругих соединений, сосредоточившись на подсистемах 
приводов: 
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В самом простом случае практический интерес представляет анализ области существова- 
ния стационарных положительных значений скорости подачи и скорости вращения шпинделя в 
зависимости от сил, приложенных со стороны приводов. Для этого уравнение статики для систе- 
мы (7) необходимо разрешить относительно координат состояния. После очевидных преобразова- 
ний и подстановок получаем уравнение для скорости главного движения: 


В 69 + (№, К, — ВМ, ) в. +, М, =0, (8) 


(9) 





Необходимым условием существования сверления является неотрицательность дискрими- 
нанта. Это позволяет определить возможные значения скорости вращения в зависимости от сил, 
прикладываемых со стороны привода: 
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Рис. 4. Необходимое условие существования сверления [по уравнению (10)] (а) 
и частота вращения привода главного движения [действительные корни уравнения (9)] (6) 


66 


Вестник ДГТУ. 2012. № 1(62), вым. 








Значение М? = 0,7353 крутящего момента на приводе главного движения соответствует 


(при данной осевой нагрузке на сверло) бифуркационному значению, при котором сверление ста- 
новится возможным. Диаграммы стационарных точек -— зависимость скорости вращения шпинделя 


от М, и РЁ, приведены на рис. 5. 
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Рис. 5. Зависимость скорости вращения шпинделя [уравнение (9)] от М, и р, 


Следующий вопрос, который должен быть исследован -— устойчивость стационарных точек. 
Линеаризуя (7) в стационарных точках для значений М’ и Е, ‚ получаем системную матрицу: 
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Для значений силовых воздействий М! =0,76 Н.ми Е, =3000 Н определяем стационар- 
ные частоты вращения ‹>’” = 0,4145 и 5? =0,7174, скорости у =0,0101 и у? =0,0102 и, та- 
ким образом, получаем явный вид системных матриц: 

и 21,9 2090 @) _ —12,0 1208 


(12) 
115 —3145 0,5 —3048 


Матрица А“ имеет собственные числа ^(' = 22,8 и /.0 =-3146 и соответствует неустой- 


чивому положению равновесия (нижняя ветвь на рис. 4,6), матрица А) имеет собственные чис- 
ла А.) =—11,8 и А) =—3049 и соответствует устойчивому положению равновесия (верхняя 


ветвь на рис. 4,6). 
Фазовая траектория нелинейной системы сверления, полученная численным моделирова- 
нием и подтверждающая приведенные выше рассуждения и выкладки, показана на рис. 6. 
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Рис. 6. Фазовый портрет системы сверления для № и =2,75 Н.м, Е =500Н 


и координатами { в, у, (рад/с, мм/с) устойчивой 43,31: 1,4} 


и неустойчивой {0,28; 0,6} точек соответственно 


Заключение. Изложена процедура построения математической модели спецстанка для сверле- 
ния глубоких отверстий в труднообрабатываемых и гетерогенных материалах. В ходе проведения 
работ были получены следующие результаты: 

— получена упрощенная математическая модель динамики сверления на спецстанке для 
сверления глубоких отверстий; 

— при анализе модели определены условия возможности сверления (существование и ве- 
личина бифуркационного параметра) в зависимости от характеристик преобразующей системы 
спецстанка и силовых воздействий приводов подачи и главного движения; 

— построены области существования сверления в зависимости от силовых воздействий 
приводов подачи и главного движения; 

— проведено имитационное моделирование и получен пример фазового портрета, под- 
тверждающего правильность проведенных аналитических вычислений. 
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ОУМАМТС МООЕЕ ОРГ ОЕЕР НОГЕ ОВКТТМС ЗРЕСТАЕ ТОО РОК НАКО МАТЕКТАТ$ 
ГАКСЕ ОТАМЕТЕК ТОВЕ РЕАТЕ ТВЕАТТМС 


Т.$. ОМОУКО 

(Роп Чае Тесптса! Упмег$Ку), 
Р.Р. УЕВЕЗНСНЕТТМ 
(Упмег$Ку ог Мгойта, Ч$А), 
А.О. СОКУАМО\У 

(Роп Чае Тесйптса! Утмег$Ку) 


Реаиге5 о! {пе татетайса!/ тодЕ/тд о! ата ю0/5 тодегпеа Гог а//!тд деер Роз аге 5 ед. А таетаг- 
Гса/ тоде! о! а!И/!тпа, диезйоп$ о! $аКу, апа ехт$епсе о! те ап/!тд ргосех$ аге сопэ!АЕгед. 
Кеуигога$: деер ап/!тд, тодеЕ/тд, дупат/с$, ат! ооЕ. 
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